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Severe declines in honey bee populations have made it imperative to understand key factors impacting 
honey bee health. Of major concern is nutrition, as malnutrition in honey bees is associated with immune 
system impairment and increased pesticide susceptibility. Beekeepers often feed high fructose corn syrup 
(HFCS) or sucrose after harvesting honey or during periods of nectar dearth. We report that, relative to 
honey, chronic feeding of either of these two alternative carbohydrate sources elicited hundreds of 
differences in gene expression in the fat body, a peripheral nutrient-sensing tissue analogous to vertebrate 
liver and adipose tissues. These expression differences included genes involved in protein metabolism and 
oxidation-reduction, including some involved in tyrosine and phenylalanine metabolism. Differences 
between HFCS and sucrose diets were much more subtle and included a few genes involved in carbohydrate 
and lipid metabolism. Our results suggest that bees receive nutritional components from honey that are not 
provided by alternative food sources widely used in apiculture. 

Honey bees are vital members of natural and agricultural ecosystems worldwide. In the United States, the 
Western honeybee (Apis mellifera) contributes more than 15 billion dollars to the agricultural industry 
annually 1 . It is therefore of serious concern that honey bee populations have declined steadily in the 
United States, with dramatic losses of colonies starting in 2006 associated with colony collapse disorder (CCD) 2 4 . 
These losses have intensified the need to understand factors that impact honey bee health. 

Central to honey bee health is nutrition 5 . Malnutrition in honey bee colonies can result from maintaining 
densities of colonies that are too high for available flora or placement of colonies for pollination of crops that are 
deficient in pollen or nectar or have low nutritive value 5-6 . Poor nutrition can make bees more susceptible to 
pesticides 7 and lead to a compromised immune system making bees more vulnerable to diseases 8 . 

The principal natural carbohydrate source of honey bees is nectar, which is collected from flowers, transported 
to the hive and converted to honey for storage. This conversion involves reducing the water content to 16-20% 
and adding glandular secretions that contain microorganisms and enzymes, including amylases, glucose oxidases 
and invertases 5,9 . These increase acidity and convert the sucrose in nectar into glucose and fructose 9 . The final 
constituents of honey vary depending on the nectar source but are mainly fructose (30-45%), glucose (24-40%) 
and sucrose (0.1-4.8%), as well as trace amounts of other disaccharides, vitamins, minerals, amino acids and a 
variety phenolic compounds 10 . 

Adult honey bees use honey as fuel for energy-intensive flights and colony thermoregulation 9 . Unlike larvae, 
adults have low levels of abdominal lipids and cannot survive for long periods of time without a carbohydrate 
source. A continuous supply of sugar is particularly important for foraging honey bees, because they have a diet 
that is mainly carbohydrate-based 11 . Compared to younger bees that specialize on performing tasks inside the 
hive, foragers also have a higher metabolic rate 12 and lose over half their abdominal lipid stores prior to starting to 
forage 13 . 

Beekeepers often provide supplemental carbohydrates in the form of high fructose corn syrup (HFCS) or 
sucrose following the harvesting of honey or during periods of nectar dearth. Supplementing with HFCS became a 
widespread practice following early studies that showed acceptable honey bee survival 14 and equivalent honey 
production and long-term productivity relative to honey feeding 15 . In addition, HFCS has a fructose-to-glucose 
ratio similar to honey, with the most common bee feed formulation composed of 55% fructose and 42% glucose 16 . 
HFCS is also less expensive than sucrose and is less labor-intensive to administer as food because it comes in liquid 
form 17 . 
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Figure 1 | Differences in gene expression in honey bee fat body caused by diets of honey, sucrose or high-fructose corn syrup (HFCS). (A) Number of 
differentially expressed genes (DEGs) for each diet comparison (FDR<0.1). (B) Number of DEGs that overlap among the diet comparisons. (C) Gene 
Ontology categories significantly enriched (P<0.004) for each diet comparison. 



However, questions regarding the suitability of HFCS for honey 
bees have arisen, in part because of CCD and because of research 
showing HFCS may have deleterious metabolic effects in mammals 18 . 
Artificial honey produced exclusively from HFCS has a different 
carbohydrate profile (contains fructosyl-fructose) than artificial 
honey produced from sucrose 19 and has been shown to decrease 
spring brood and wax production relative to sucrose 20 . Moreover, 
there is growing evidence that constituents in natural honey, absent 
from sucrose and HFCS, positively affect the honey bee's xenobiotic 
detoxification system 21,22 . These results suggest that honey, sucrose 
and HFCS may impact honey bee physiology and health differently. 

We explored this issue further with whole-genome transcrip- 
tomics to comprehensively survey the effects of honey, sucrose and 
HFCS on fat body gene expression. The fat body is a multifunctional 
organ responsible for nutrient storage, energy mobilization and the 
production of antimicrobial peptides. Nutrient storage and mobil- 
ization are coupled to hormonal signals that include insulin and 
adipokinetic hormone to fulfill ongoing physiological demands 23 . 
In adult honey bees, the fat body is known to be transcriptionally 
responsive to nutritional manipulations and manipulations that 
affect aging and health 24,25 . We focused on the fat body to study 
the effects of different dietary carbohydrate sources on the expression 
of genes involved in hormonal signaling, nutrient storage, energy 
metabolism, and immune function. In this study, we used older bees 
(18-21 days old) because their diet is primarily carbohydrate-based 11 
and because older bees have been shown to be the primary consu- 
mers of carbohydrate supplements inside the hive 5 . 

Results 

Measurements taken daily throughout the week-long trials showed 
similar levels of food consumption for bees fed honey, HFCS or 
sucrose (0.040 ± 0.001, 0.036 ± 0.002, 0.032 ± 0.003g/bee/day, 
respectively, F=4.26 P=0.055). Mortality also did not vary between 
diets (F=0.57 P=0.59) and was between 0-7% across all cages. 

RNA-sequencing (RNA-seq) was performed to examine the effect 
of each diet treatment on fat body gene expression. In total, 5 sample 
pools were sequenced per diet treatment and colony replicate 



(N=10). Initial examination of our results revealed that one of our 
colony replicates was heavily infected with deformed wing virus 
(DWV). For Colony A, an average of 37.06 ± 8.66% reads aligned 
to the DWV genome sequence compared to only 1.29 ± 2.77% for 
Colony B. By contrast, an average of 53.81 ± 0.08% of the RNA-seq 
reads mapped to the honey bee genome for Colony A, while for 
Colony B the average was 87.72 ± 0.03%. Due to this difference, 
we analyzed each colony separately. We explored the effects of diet 
treatment using multidimensional scaling (MDS) plots (using the log 
fold change of the top 100 genes) generated for each colony. These 
analyses indicated that chronic feeding of sucrose and HFCS elicited 
transcriptomic profiles distinct from bees fed a honey diet (Figure 
S1A and SIB). This pattern was observed for both colonies, indi- 
cating similar transcriptomic responses to diet treatments occurred 
regardless of differences in apparent viral load. 

To probe more extensively for diet effects on gene expression, we 
analyzed the results from both colonies together, assessing the main 
effect of diet treatment with colony as a blocking factor. Honey 
elicited hundreds of differences in gene expression relative to 
HFCS and sucrose (Figure 1A). There were 104 genes differentially 
expressed (FDR <0.1) in bees fed honey or HFCS and 220 genes 
differentially expressed between bees fed honey or sucrose. By con- 
trast, differences between HFCS and sucrose diets were much more 
limited with a total of 8 genes differentially expressed between these 
two groups. Gene-wise comparisons show a substantial overlap (64 
genes) between the honey vs. sucrose (29.1% overlapped) and honey 
vs. HFCS (61.5% genes overlapped) indicating strong similarities 
across these gene lists (Figure IB). 

Class prediction analyses using the support vector machine algo- 
rithm 26 revealed that diet-induced fat body gene expression changes 
were robust and consistent across samples. Class membership was 
predicted correctly with 96% (Honey vs. Sucrose), 97.5% (Honey vs. 
HFCS) and 100% (HFCS vs. Sucrose) accuracy, corresponding to 
sensitivity values (true positives identified) between 0.92-1 and spe- 
cificity values (true negatives identified) of 0.95-1. Top predictors for 
each diet comparison are shown in Figure S2. Notable among top 
predictor genes for the honey-based comparisons were glutathione S 
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transferase 03 (GB44803) and pale (GB40967), which are associated 
with xenobiotic detoxification and tyrosine metabolism, respectively. 
Maltase Bl (GB54549) and two other genes involved in energy meta- 
bolism, GB50596 (GO term: oxido-reductase activity) and GB48029 
(GO term: acyl carnitine transporter activity), were top predictors for 
the HFCS vs. Sucrose comparison. 

Consistent with class prediction analyses, Gene Ontology (GO) 
enrichment analyses showed that the lists of genes upregulated by 
honey were enriched for genes involved in amino acid metabolism 
and oxidation reduction (Figure 1C, Table SI), especially phenylala- 
nine and tyrosine metabolism. These included pale, henna 
(GB48022) and homogentisate 1,2-dioxygenase (GB53288). Relative 
to sucrose, honey also upregulated the gene flavin monooxygenase 1 
(GB42239), which was associated with oxidation reduction and 
alkaloid detoxification 27 . By contrast, sucrose upregulated genes that 
were associated with axonogenesis, anion transport, and several tran- 
scription factors associated with organ development, such as double- 
sex (GB55036), knot (GB42304) and tolkin (GB52106), while HFCS 
upregulated the transmembrane receptors, domeless (GB42244) and 
tyramine receptor (GB47385) (Figure 1C, Table SI). No GO terms 
were enriched for the 8 genes differentially expressed between HFCS- 
and sucrose-fed bees. 

To gain further insights into the biological significance of these 
diet-induced differences in fat body gene expression, we compared 
our results to three previously published transcriptomic (microar- 
ray) experiments on nutritional aspects of behavioral maturation in 
honey bees 24 . Behavioral maturation in honey bees involves a switch 
from a high protein to a high carbohydrate diet and a loss of approxi- 
mately 50% of fat body lipids 13 prior to the shift from working in the 
hive to foraging. This behavioral and physiological shift also includes 
a reduction in the blood titers of the lipid storage protein vitellogenin 
(Vg) 28 which has been associated with the promotion of honey bee 
longevity 29 and immunity 30 . We compared our results to three micro- 
array experiments: 1) Maturation (hive bees compared to foragers); 
2) vg knockdown (vgRNAi compared to control); and 3) Diet [bees 
fed a high protein diet (45% pollen, 45% honey, 10% water) com- 
pared to sucrose (50% w/v)]. We detected significant overlap 
between previously reported changes in fat body gene expression 
that occur during the hive to foraging transition 24 and our present 
Honey vs. Sucrose and Honey vs. HFCS gene lists (Table 1). 
Surprisingly, the fat body gene expression profile of the nutritionally 
enriched hive bees was more similar to that of the sucrose-fed bees, 
while the profile of the more nutritionally deprived foragers was 
more similar to that of the honey-fed bees (Table 1). GO enrichment 
analyses showed that the subset of genes overlapping the Honey vs. 
Sucrose or Honey vs. HFCS diets and the maturation-related experi- 
ment were associated with protein metabolism and oxidation reduc- 
tion (Table S2). There also was a significant overlap between the 
Honey vs. Sucrose and Honey vs. HFCS gene lists and the gene list 
from the vg RNAi experiment 24 ; however, the log fold changes for 
overlapping genes were not significantly correlated indicating gene 
lists were not directionally similar (Table 1). There was no significant 



overlap between our Honey vs. Sucrose and Honey vs. HFCS gene 
lists and the Diet experiment 24 , but the genes that did overlap showed 
a significant positive correlation that suggested honey's effects were 
directionally concordant to those of the Pollen + Honey treatment 
(Table 1). 

Discussion 

A honey diet elicited a transcriptional profile distinct from sucrose 
and HFCS diets. These differences were present in two different 
honey bee colonies, with vastly different viral loads, indicating the 
impact of honey on fat body gene expression is robust. These results 
suggest that constituents in honey differentially regulate physio- 
logical processes and that sucrose and HFCS may not be equivalent 
nutritional substitutes to honey. 

Gene Ontology enrichment analyses showed honey upregulates 
genes associated with processes such as "aromatic amino acid family 
metabolic process," as well as "oxidation reduction." Among the 
genes in these categories were orthologs for the Drosophila melano- 
gaster genes pale and henna, which are related to phenylalanine and 
tyrosine metabolism. These amino acids have been linked to the 
production of neurotransmitters 31 , and in the case oipale to immune 
responses to infection 32 . Honey additionally upregulated the gene 
glutathione S transferase 03, whose activity is known to be induced 
by plant compounds and to have toxicological significance in the 
presence of pesticides 33 . HFCS and sucrose relative to honey resulted 
in the upregulation of different biological processes. Sucrose, for 
example, upregulated processes such as axonogenesis but it is 
unlikely that axonogenesis is upregulated in our fat body samples; 
rather this GO category reflects upregulation of signaling pathways 
that play different roles in different tissues. HFCS upregulated the 
transmembrane receptors domeless and tyramine receptor suggesting 
differences in JAK-STAT signaling and tyrosine signaling between 
HFCS and honey. 

Sucrose and HFCS elicited a remarkably similar fat body tran- 
scriptional response. This result is consistent with previous studies 
showing no differences in colony productivity due to these diets 15 but 
contrasts with findings showing differences in wax production and 
honey bee survival due to sucrose or HFCS 14,20 . Our results suggest 
that older bees may not be sensitive to increased fructose consump- 
tion because we detected little evidence that reflects the types of 
changes in hormonal signaling, energy metabolism and nutrient 
storage associated with high fructose corn syrup and increased fruct- 
ose consumption in mammals 18,34 . Future research should test 
whether greater differences in gene expression due to sucrose or 
HFCS feeding are observed in young or middle-aged bees, which 
perform brood rearing and comb building tasks respectively, or 
whether the few differences in gene expression observed in this 
experiment can account for differences in colony performance in 
more natural conditions. 

To understand whether the gene expression differences associated 
with a honey diet relate to maturation-related physiological changes, 
we compared our results with previously published fat body micro- 



Table 1 | Enrichment and correlation analyses for honey-related contrasts and previously published maturation-related microarray experi- 
ments. Pairwise comparisons are for genes differentially expressed at an FDR<0. 1 . The number of genes compared is shown next to each 
experiment name in parentheses. RF is an enrichment factor for the number of genes that overlapped and P values indicate significant 
enrichment (hypergeometric test). Correlation results were calculated with genes that overlapped between pairwise comparisons 

Nurse vs. Forager (21 54) Control vs. vg RNAi (2823) Pollen + honey vs. sucrose (2449) 
Experiment 

Overlap RF Correlation (r) Overlap RF Correlation (r) Overlap RF Correlation (r) 

Honey vs. Sucrose 86 1.306 -0.261 (P=0.03) 162 1.906 0.051 (n.s.) 70 0.935 (n.s.) 0.662 

(220) (P=0.002) (P<0.0001) (P<0.0001) 

Honey vs. HFCS 45 1.446 -0.1 27 (n.s.) 75 1.866 -0.1 29 (n.s.) 40 1.130 (n.s.) 0.645 

(104) (P=0.002) (P<0.0001) (P<0.0001) 
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array studies 24 . These comparisons showed significant enrichment 
between honey-based gene lists and genes differentially expressed 
between hive bees vs. foragers, as well as with gene expression differ- 
ences associated with vitellogenin RNAi treatment. Functional ana- 
lyses suggest that shared changes in gene expression were related to 
protein metabolism and oxidation reduction, suggesting these pro- 
cesses are responsive to direct diet manipulations and maturational 
changes. Contrary to our expectations, we found the transcriptional 
profile of honey- fed-bees resembled the less well-nourished foragers 
rather than the more well-nourished nurse bees. These results sug- 
gest that there may be compounds in honey that modulate honey bee 
physiology towards a forager-like state. 

We did not find significant overlap between our Honey vs. Sugar 
differentially expressed gene list and the previously published list 
from bees that received either a diet of Pollen + Honey or sucrose. 
This result is surprising because the former contrast was embedded 
within this diet-related microarray experiment. This indicates that 
pollen is largely responsible for the gene expression changes in the 
diet microarray experiment and that those changes are separate from 
those elicited by honey in our experiment. Lack of statistical enrich- 
ment with pollen- induced changes partially reflects the relatively 
low level of protein in honey. In addition to the use of different 
experimental platforms (RNA-seq vs. microarray), there were also 
differences in the age of the bees assayed in each experiment: the diet- 
related microarray experiment investigated the effect of diet treat- 
ment on younger bees able to digest pollen, while our experiment 
assayed older bees with a decreased capability to digest pollen. Thus, 
the transcriptional differences elicited by a honey diet cannot be 
directly attributed to pollen traces in honey. 

Our goal was to perform a broad unbiased survey for the effects of 
honey, sucrose and HFCS on honey bee physiology. Our result that 
honey - but not sucrose or HFCS - upregulates genes associated with 
protein metabolism and oxidation reduction is indicative that honey 
elicits health-related physiological differences. We performed our 
experiment using older bees that typically consume sugar solutions 
inside the hive and do not digest pollen; therefore, constituents in 
honey may provide critical nutritional components and inducers that 
are otherwise absent in this age group. Previous research has already 
identified honey constituents that upregulate detoxification path- 
ways in the gut 22 ; our results further show honey induces gene 
expression changes on a more global scale. These changes may have 
toxicological relevance under natural conditions in contemporary 
agroecosystems, where bees are routinely exposed to toxins and 
pesticides. 

Methods 

Bees. We used bees from honey bee colonies from the University of Illinois Bee 
Research Facility, Urbana, IL, maintained according to standard beekeeping 
practices. The bees were a mixture of European subspecies typical of this region. To 
minimize genetic variation within a replicate, we used adult worker bees from a 
colony derived from a queen inseminated by single male; due to haplodiploidy, these 
bees were related to each other by an average coefficient of relatedness of 0.75. The 
experiment was replicated in two independent trials, each time using bees from 
different, unrelated, colony. 

Feeding Trials. We used adult bees between the ages of 18-21 days old. These are 
older bees that readily consume various carbohydrate sources in the hive 5 . To obtain 
focal bees we removed honeycomb frames containing pupae, placed them in an 
incubator (34°C/30% RH), marked newly emerged one-day-old bees with a spot of 
paint (Tester's Paint, Rockford, IL, USA) on the dorsal surface of the thorax and 
reintroduced the marked bees into their natal colony; this was repeated three 
consecutive days to obtain a base population of >500 marked bees in the hive. Focal 
bees were collected when they were 18-21 days old, placed into Plexiglas cages (10 X 
10 X 7 cm; 15 bees per cage) and assigned a diet treatment in the laboratory. Diet 
treatments consisted of 50% (w/v) honey, 50% high fructose corn syrup {HFCS 55) or 
50% sucrose ad libitum for each cage replicate (N = 3 cage replicates per treatment). 
All cages were kept in a 29°C incubator. Consumption and mortality were monitored 
daily for 7 days. After this time period, bees were flash-frozen and stored at — 80°Cfor 
analysis. 



Dissections and RNA extractions. Dissections were performed by incubating each 
abdomen in chilled RNA-later ICE (Ambion) at — 20°C for a minimum of 16 h. The 
gut and ventral tissue were then removed and RNA extraction performed on the fat 
body and adhering dorsal cuticle {RNeasy kit, Qiagen, with a DNAse treatment). 

cDNA library construction and RNA-sequencing. We constructed cDNA libraries 
using pooled total RNA from fat body tissue from 3 individual bees. Pooling was 
performed to minimize sample variability within treatment groups. Each cDNA 
library was prepared with 1.5 ug of pooled total RNA and constructed using the 
NEXTflex™ Directional RNA-seq kit (Bioo Scientific) with an added mRNA 
purification step using Dynabeads® Oligo(dt) 25 (Invitrogen). Library concentrations 
were quantified using Qubit® fluorometric quantitation and by quantitative real-time 
PCR using Kapa Library Quant kits (KapaBiosys terns). Average fragment size and 
overall quality was evaluated with the Agilent 2100 Bioanalyzer platform and an 
Agilent High Sensitivity DNA kit. For sequencing, all libraries were diluted to a 6 nM 
concentration. In total, we sequenced 5 libraries per treatment per colony, for a total 
of 30 libraries. Ten libraries were sequenced per lane (2-4 libraries per diet treatment/ 
per lane) with an Illumina HiSeq2000 instrument and were sequenced as single-end, 
100 nt reads. 

Bioinformatics. RNA sequencing generated an average of 19,200,920 reads per 
library. Reads from each library were aligned against the Apis mellifera genome, 
Assembly 4.5 35 using Tophat 36 and read counts were generated for genes using HTSeq 
v0.5pv2 (http://www-huber.embl.de/users/anders/HTSeq) and the Apis mellifera 
Official Gene Set, version 3.2. Reads that did not map uniquely or that mapped to 
genomic locations outside genes were not included in analyses for differential 
expression. We identified differentially expressed genes (DEGs) using a generalized 
linear model in the EdgeR package, version 3.2.4 37 . 

Because the number of differentially expressed genes in this study was in the 
hundreds instead of the thousands, we performed an additional analysis to test the 
robustness of the results. We compared differences in fat body gene expression 
between two subsets of the entire group of individuals fed on honey (n— 5 per subset), 
permuting the samples allocated to each subset a total of 100 times. We found only a 
few genes [2.07 ± 0.344 genes (FDR<0.1)] to be differentially expressed in these 
comparisons, suggesting that the differences we report between diet treatments are 
reliable. 

Identification of Deformed Wing Virus Infection. Reads that did not align to the 
Apis mellifera genome were queried using BLAST on the NCBI website. BLAST 
results identified unaligned reads from Colony A as deformed wing virus (DWV) 
sequences. We further validated DWV infection by aligning all reads from each 
sample to the DWV genome (ftp.ncbi.nlm.nih.gov/genomes/Viruses/Deformed_ 
wing_virus_uidl4891), using Bowtie 2. 

Class Prediction Analyses. To identify the most robust and consistent changes in fat 
body gene expression caused by each diet, we performed class prediction analyses 
using the support vector machine algorithm with 5-fold cross-validation. Class 
prediction analyses were performed using CMA 26 , with normalized read counts for 
genes with a P value <0.05. Top predictors for each diet contrast were selected based 
on their ranking as a top classifier in at least 10 of the 50 iterations. 

Gene Ontology Analysis. Inferences for major functional themes for each DEG list 
were drawn from GO enrichment analyses using the DAVID Bioinformatic 
Resources 6.7 functional annotation tool 38 . These analyses were performed using the 
Drosophila melanogaster orthologs associated with each DEG list. Statistical analyses 
of enrichment were performed with a hypergeometric test, using the number of honey 
bee genes with annotated Drosophila orthologs as the reference. 

Comparisons to Previous Microarray Studies. Ament et al. (2011) previously 
published the results of microarray experiments investigating fat body gene 
expression related to nutritional aspects of behavioral maturation. We compared our 
results to three of them: 1) Maturation (hive bees compared to foragers; 2) vitellogenin 
knockdown {vitellogenin dsRNA compared to control; and 3) Diet: [bees fed a high- 
protein diet (45% pollen, 45% honey, 10% water) compared to sucrose (50% w/v)]. 

To determine whether the lists of DEGs contained significant levels of overlap, we 
calculated an enrichment factor (RF) by dividing the observed number of overlapping 
genes by the expected number. The expected number of overlapping genes was 
calculated by multiplying the length of each DEG list and then dividing this value by 
the total number of genes included in the analyses 39 . An RF value greater than 1 
indicates the observed number of overlapping genes is greater than the expected 
value, thus demonstrating enrichment. Significant enrichment was determined using 
a hypergeotmetric test (1 -tailed) with the p-hyper function in R. To determine 
whether the overlapping genes between two lists were directionally concordant in the 
predicted direction, we compared the log-fold changes for each gene list and assessed 
significance with a Pearson's correlation. A positive and significant Pearson's cor- 
relation shows the majority of overlapping genes between two gene lists show con- 
cordant gene expression changes. 
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